Reexamination for the calculation of elliptic flow and other fourier harmonics 
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We have argued that the azimuthal symmetry and asymmetry components in fourier expansion 
of particle momentum azimuthal distribution, v n (n=0, 1, 2, ...), should be calculated as an average 
of cos(n(j>) first over particles in an event and then over events (event-wise average) rather than 
"an average over all particles in all events" (particle- wise average). In case of large centrality 
(multiplicity) bin the particle-wise average is not accurate because the influence (fluctuation) of 
particle multiplicity was not taken into account. 
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Elliptic flow (proportional to v%) and other harmonics 
(proportional to v n , n=0, 1, 3, 4, ...) as fourier expansion 
coefficients of azimuthal distribution of particle momen- 
tum are highly sensitive to the eccentricity of the early 
created fireball in ultra-relativistic heavy ion collisions. 
The expected phase transition to Quark-Gluon-Plasma 
(QGP) should have a dramatic effect on those harmonics. 
The consistency between experimental data of V2(pt) at 
mid-rapidity and the corresponding hydrodynamic pre- 
dictions are regarded as an evidence of the production of 
partonic matter in ultra-relativistic nucleus-nucleus col- 
lisions [lj . The elliptic flow of high px particles may be 
related to the jet fragmentation and parton energy loss 
[2|, which are usually not included in the hydrodynamic 
calculations. This hydrodynamic model J3|] overestimates 
V2(pt) in the pr > 1.5 GeV/c region [4J. That is re- 
garded, together with the discovery of jet quenching Q, 
as a strong evidence of sQGP formation in relativistic 
nucleus- nucleus collisions at RHIC. 

The elliptic flow and other harmonics have become an 
interesting topic in the field of ultra-relativistic nucleus- 
nucleus collisions. A lot of experimental data from RHIC 
have been published [fj, 0, Q . Consequently the micro- 
scopic transport model studies are also widely progress- 
ing [3, [H], [13, E2, O as well as the abundant hydrody- 
namic investigations. 

Refs. [l4J 1 151 are two well known pioneering papers in 
this field. In (l5| | the investigation was starting from the 
triple differential distribution of 
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[1+ J2 2u„cos[ra(<£-# r )]]. (1) 



Then it was defined that: "... v n — (cos[n(<f> — ^V)]), 
where () indicates an average over all particles in all 
events. For the particle number distribution, the co- 
efficient vi is (p x /pt) and v 2 is ((p x /p t ) 2 - (p y /p t ) 2 )." 
(This kind of average will be indicated as particle-wise 
average hereafter.) Later that is widely accepted in the- 
oretical calculations and becomes a common convention: 



1 and n-th harmonic is 
3, ... ), here () indicates 



the zeroth harmonic is i>n 
v n = {cos[n((j> - \]> r )]) (n=l 
particle- wise average. 

As mentioned above that the elliptic flow and other 
harmonics are important and widely studied for over a 
decade. But we have to review the calculation of v n in 
this paper: does the v n should be calculated as an average 
of cos(n(j>) first over particles in an event and then over 
events (event-wise average) or calculated as an average 
of cos(n(f)) over all particles in all events (particle- wise 
average). To this end, we first rederive the elliptic flow 
and other harmonics starting from the particle number 
(multiplicity) distribution. For that deduction one has 
first to define the reaction plane. 

In theory, if the beam direction and impact parameter 
vector are fixed, respectively, at the p z and p x axes, then 
the reaction plane is just the p x —p z plane jl4| . Therefore 
the reaction plane angle (\& r ) between reaction plane and 
the p x axis [14J introduced for the extraction of elliptic 
flow in the experiments [15J ] is zero. The particle az- 
imuthal angle ((/>) in momentum space is measured with 
respect to the reaction plane, which is consistent with 
the definition in [l5|]. This particle azimuthal angle (4>) 
is just the angle spanned by pr relative to the p x axis as 
shown in Fig. [T] 




FIG. 1: The particle momentum p, in the Cartesian coordi- 
nate system, \p x ,p y ,Pz]- 
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However, in experiment, the reaction plane is different 
event by event, thus the experimental measurement of 
elliptic flow is not trivial. One has to invoke a complex 



reaction plane identification method (two-particle corre- 
lation method) Hal , cumulant method [16J], or Lee- Yang 
zeroes method [17| . In all of these methods a quantity has 
to be first constructed event by event: the event plane 
in [15| , cumulant expansion of weighted n-th transverse 
event- flow vector in [16j |. and a generating function in 
[17| . Then a corresponding average over measured events 
has to be taken. Therefore the experimental extraction 
of elliptic flow parameter is event-wise average indeed. 

The particle number (multiplicity) distribution can be 
expressed as 
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in the cylindrical coordinate system after substituting p z 
by y (rapidity) and using the relation of dy/dp z — 1/E 
[18j |. Then the normalized particle multiplicity distribu- 
tion is 



1 



d 3 N 



N ev p T dydp T c 



N ev = dy p T dp T 



d 3 N 



prdydpTC 



(3) 



where N ev is the event multiplicity and the integrals are 
taken over entire range of the variables. 

This normalized particle multiplicity distribution in 
Eq.Q is a three dimensional (y,pT,4>) distribution func- 
tion. The dimension can be reduced by integrating over 
a certain variable [lj|. For example, to study the v n as 
a function of rapidity y, v n {y), one should take integral 
over px, then above three dimensional distribution func- 
tion reduces to a two dimensional distribution function 
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In numerical calculations a small rapidity interval, Ay, 
is used instead of single y value, the corresponding nor- 
malized particle multiplicity distribution becomes 



— [dy' f Pl d Pl d3N 
N evp J p J prdy'dprc 
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d(j) N evp J p J p T dy'dp T 



where 



N eV p = dy' prdpT / i 
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and 
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Note that, the normalization factor N evp here is differ- 
ent from N ev in Eq. ([3]) in the range of integral over 
y and is denoted as constrained event multiplicity. In 
the above equation f((j>) is the normalized distribution 
density function of <f> and the N evp f((j)) is the number of 
particles emitted into d<f> at azimuthal angle cf> without 
constraint on px but y is constrained in Ay. This can 
be constructed experimentally and/or calculated theoret- 
ically. Here we have taken v n (y) as an example but it is 
similar for the v n (px)- 

Since f((f>) is 2ir periodic and even function of </>, it can 
be expanded by a fourier series [2CJ as 
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cos(n4>) (n = 1,2, ...). 



(7) 



In above equation, cos(n<j)) denotes the average of 
cos(n(f>) over particles in a single event and can be calcu- 
lated both experimentally and theoretically. The first 
expansion term, in Eq.®, is a circle (isotropic), sec- 
ond term (cos(4>)) a leave (directed flow), third term 
(cos(2(j))) a four-leaved rose (elliptic flow), and the (n+1)- 
th term (cos(n(j))) a multi-leaved rose (n-th harmonic) 

As the event multiplicity fluctuates event-by-event, one 
always has to generate multiple events and to take aver- 
age over events generated. So the v n in Eq.([7]) should 
be 



v e n = (cos(n(j)))ev (n = l,2, ...), 



(8) 



where (...) e v means an average over all events. A super- 
script , "e" , is here added on v n to identify this average 
as the event-wise average. Applying recursion formula of 
cosine function 20] 



cos(n<p) — 2cos((n — l)<p)cos<p — cos((n — 2)(f>) (9) 
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to Eq. (|5J|, v^ can be expressed as 
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Above deductions argued that the v n is an average of 
cos(n<fi) first over particles in an event and then average 
over events (event-wise average) rather than a particle- 
wise average ( "an average over all particles in all events" 



fig]) of cos{n<p). The particle- wise average of cos(nq 
can be expressed as 



(cos{n(/))N evp } ev /(N evp } 



v evp/ev 
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(see Appendix). This is obviously different from the 
event- wise average. Only if the cos(n(j)) is independent 
of constrained event multiplicity N evp (i. e. constrained 
multiplicity plays no role in the average) the v? reduces 
to v!^. In fact, the cos(n<fi) and N evp correlate (even 
negatively correlate) with each other. This is because 
the larger event multiplicity is due to more central col- 
lision (larger overlap zone between colliding nuclei) and 
the larger overlap zone in turn leads to less azimuthal 
asymmetry. The particle- wise average does not take the 
influence of event multiplicity into account, thus it is in- 
accurate from physics point of view. Of course, for very 
narrow multiplicity (centrality) bins the particle- wise av- 
erage is not very problematic. However for the wide mul- 
tiplicity bins (such as the 0-40% most central Au+Au col- 
lisions studied in [2l| and /or the minimum bias Au+Au 
collision data given in [22|]) that correction is not negli- 
gible. 

We have applied the parton and hadron cascade model, 
PACIAE [23| , to calculate the charged hadron vi (77) in 0- 
40% most central Au + Au collisions at • v /siviv=200 GeV 
by the method of event-wise and particle-wise averages 
separately. The v%{r\) calculated by event-wise average 
is about 10% larger than v\(r\) calculated by particles- 
wise average. This means that {cos(n<f))N evp ) ev is smaller 
than {cos{n4>)) ev {N evp ) ev , i. e. there is a negative corre- 
lation between cos(n(j>) and N evp . This is consistent with 
the physical picture given in the last paragraph. 

In summary we have rederived the elliptic flow and 
other harmonics in fourier expansion of the particle mo- 
mentum azimuthal distribution starting from the parti- 
cle number (multiplicity) distribution in the momentum 
space. It is argued that the parameter oin—th harmonic, 
v n , should be calculated as an average of con(n<p) first 
over particles in an event and then over events (event- 
wise average). This is different from the particle- wise 
average ( "an average over all particles in all events" ) of 
cos(n(f>). In case of large centrality (multiplicity) bin the 



particle-wise average is not accurate because the influ- 
ence (fluctuation) of particle multiplicity is not taken into 
account. 

Appendix 

In this appendix the variables are assumed to be discrete 
for simplicity. Then the particle- wise average of cos(ncj)) 
can be expressed as 



J2[cos(n(j)))k 

_k 

k 



(13) 



where J2 denotes the sum over all particles in all events. 

k 

If the summation, both in numerator and denominator, 
groups according to event (denoted by j) then 
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where E is the sum over events and ^2 denotes the sum 

3 *ej 

over particles in the j-th event. According to the defini- 
tion of average over particle in a single event [20| the v? 
can be further expressed as 



J2(cos(n4>))j(N evp )j 
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